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Background of Invention 

[0001] The invention relates generally to the field of well logging. More particularly, the 
invention relates to improved techniques for well logging using nuclear magnetic 
resonance and methods for analyzing the nuclear magnetic measurements. 

[0002] Oil well logging tools include nuclear magnetic resonance (NMR) instruments. NMR 
instruments can provide a wealth of information for formation evaluation that is not 
obtainable from other well logging measurements. Information provided by NMR 
measurements include the fractional volume of pore space, the fractional volume of 
mobile fluid filling the pore space, and the porosity of earth formations. General 
background of NMR well logging is described in U.S. Patent No. 6,1 40,81 7A1 , 
assigned to the assignee hereof. 

[0003] 

The signals measured by nuclear magnetic resonance (NMR) logging tools typically 
arise from the selected nuclei present in the probed volume. Because hydrogen nuclei 
are the most abundant and easily detectable, most NMR logging tools are tuned to 
detect hydrogen resonance signals (form either water or hydrocarbons). These 
hydrogen nuclei have different dynamic properties (e.g., diffusion rate and rotation 
rate) that are dependent on their environments. The different dynamic properties of 
these nuclei manifest themselves in different nuclear spin relaxation times (i.e., spin- 
lattice relaxation time ( 7" ) and spin-spin relaxation time ( T ^ )) and diffusion 
constants. For example, hydrogen nuclei in viscous oils have relatively short relaxation 
times and low diffusivity, whereas hydrogen nuclei in light oils possess relatively long 
relaxation times and high diffusivity. Furthermore, the hydrogen nuclei in free water 
typically have longer relaxation times than those in bound water. Consequently, these 
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differing NMR relaxation times can provide information on properties of the earth 
formations. 

[0004] Most NMR logging tools measure the spin-spin relaxation times ( T ^ ) to derive 
the properties of the earth formations. T ^ relaxation is often measured from a train 
of spin-echoes that are generated with a series of pulses such as the Carr-Purcell- 
Meiboom-Gill (CPMC) pulse sequence or some variants of this. The CPMC pulse 
sequence is well known in the art. See Meiboom, S., Gill, D., 1 958, "Modified Spin Echo 
Method for Measuring Nuclear Relaxation Times, "Review of Scientific Instruments, 29, 
688-91. 

[0005] Although T ^ measurements provide useful information for formation 
characterization, T information alone may not be sufficient for distinct 
characterization of earth formations, especially when different components in the 
formations have similar or overlapping T values. In order to differentiate different 
fluids (e.g., hydrocarbons versus connate water) with similar or overlapping T ' 
distributions, several differential methods have been proposed, such as the 
differential spectrum method (DSM) and time domain analysis (TDA). These methods 
takes advantages of different longitudinal relaxation times of different fluids. Two 
sets of measurements are made with different wait times (times for the spin to be 
polarized by 0 the static magnetic field). One wait time is selected that one type of fluid 
(e.g., brine that has a longer longitudinal relaxation time T ' ) would not fully relax. 
As a result, the signal magnitudes from the long T fluid would be substantially 
reduced. Subtraction of these two sets of measurements would then produce a 
difference measurement which is comprised mostly of signals from the fluids with 
long T times. However, success of these approaches relies on the selection of proper 
wait times, which requires prior knowledge of the NMR properties of the fluids in the 
formation. 

[0006] 

Recently, a magnetic resonance fluid characterization (MRF) method has been 
shown to provide more useful information. For a detailed discussion of the MRF 
method, see U.S. Patent No. 6,229,308 Bl issued to Freedman and assigned to the 
assignee of the present invention. This patent is hereby incorporated by reference. 
When T distributions overlap, the MRF method distinguishes oil and water in porous 
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media based on different molecular diffusion. For the same T ' , oil and water have 

different diffusion constants. Therefore, contributions of different fluids to the 

measured T distributions can be separated by combining measurements with 
2 

different sensitivity to diffusion. However, in order to use MRF, diffusion constants 

must be measurable. If T ^ becomes comparable with or shorter than the decay due to 

diffusion, then the diffusion constant can no longer be determined, and the MRF 

method becomes impracticable. In other words, the MRF method is applicable only 

when T values are longer than a certain critical value. 
2 

[0007] While the MRF analysis has proved to be a powerful approach, it is desirable to 
have methods that can be used to analyze fluids with not only long T or T ^ , but 
also short T or T , such as viscous fluids. Furthermore, it is desirable to have 
methods that can be used in logging tools with a low magnetic field gradient or a 
saddle point in the magnetic field. 

Summary of Invention 

[0008] One aspect of the invention relates to methods for determining properties of a 
mixture of fluids using NMRdata that include longitudinal and transverse relaxation 
information. According to embodiments of the invention, a method for determining 
properties of a mixture of fluids includes: (a) acquiring a plurality of nuclear magnetic 
resonance measurements from the mixture of fluids, each of the plurality of nuclear 
magnetic resonance measurements having a different value in an acquisition 
parameter for which at least one relaxation selected from the group consisting of 
longitudinal relaxation and transverse relaxation affects magnitudes of the nuclear 
magnetic resonance measurements; (b) generating a model of the mixture of fluids; 
(c) calculating a synthesized nuclear magnetic data set based on the model; (d) 
comparing the synthesized nuclear magnetic data set with the nuclear magnetic 
resonance measurements; and (e) adjusting the model and repeating (c) and (d), if 
difference between the synthesized nuclear magnetic data set and the nuclear 
magnetic measurements is greater than a minimum. 



[0009] 



Another aspect of the invention relates to methods for logging an earth formation 
surrounding a wellbore. According to embodiments of the invention, a method for 
logging an earth formation surrounding a wellbore includes: (a) lowering a nuclear 
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magnetic resonance instrument into the wellbore; (b) inducing a static magnetic field 
in a region of investigation; (c) generating a series of radio frequency magnetic field 
pulses in the region of investigation, and receiving signals comprising a train of 
nuclear magnetic resonance spin echoes in response to the series of radio frequency 
magnetic field pulses, wherein the generating and the receiving are repeated a 
plurality of times each with a different value in an acquisition parameter for which at 
least one of the longitudinal relaxation and transverse relaxation affects magnitudes 
of the signals; (d) generating a formation model that includes at least one component 
for a connate water phase and at least one component for an oil phase; (e) calculating 
a synthesized nuclear magnetic data set based on the formation model; (f) comparing 
the synthesized nuclear magnetic data set with the nuclear magnetic resonance 
measurements; and (h) adjusting the formation model and repeating (c) and (d), if 
difference between the synthesized nuclear magnetic data set and the nuclear 
magnetic measurements is greater than a minimum. 

[0010] Other aspects and advantages of the invention will be apparent from the following 
description and the appended claims. 

Brief Description of Drawings 

[001 1] FIG. 1 is a diagram of a typical setup for well logging. 
[001 2] FIG. 2 shows a diagram of a prior art NMR logging tool. 

[001 3] FIG. 3 is a block diagram of circuitry for producing the RF pulses and receiving the 
spin echoes. 

[0014] FIGs. 4A and 4B show the exponential recovery of signals after an inversion pulse 
or a saturation pulse. 

[001 5] FIGs. 5A and 5B show the pulse sequences for IR-CPMG and SR-CPMG, 
respectively. 

[001 6] FIG. 6 shows a flow chart of a method according to one embodiment of the 
invention. 

[0017] 

FIGs. 7 A and 7B show the results of computed porosity and T , respectively, 
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versus connate water T ^ /T values. 

[001 8] FIG. 8 shows a graph of quality of model fit to experimental data as a function of 
T /T ratios. 

Detailed Description 

[0019] FIG. 1 shows a schematic of a nuclear magnetic resonance (NMR) logging tool 30 
for investigating earth formations 31 traversed by a borehole 32 . The NMR logging 
tool 30 is suspended in the borehole 32 on an armored cable 33 , the length of which 
substantially determines the relative depth of the logging tool 30 . The cable length is 
controlled by suitable means at the surface such as a drum and winch mechanism 8 . 
Surface equipment 7 may include a processor subsystem to communicate with 
downhole equipment including NMR logging tool 30 . 

[0020] The NMR logging tool J0can be any suitable nuclear magnetic resonance logging 
device. It may be one adapted for in wireline logging applications as shown in FIG. 1 , 
or one that can be used in logging while drilling (LWD) applications. The NMR logging 
tool Jfltypically includes a means for producing a static magnetic field in the 
formations, at least one radio frequency (RF) antenna, and means for producing pulses 
of RF power to induce RF magnetic fields in the formations and for receiving the spin 
echoes from the formations. The means for producing a static magnetic field may 
comprise a permanent magnet or magnet array, and the RF antenna may comprise one 
or more RF antennas, which may be solenoid antennas, loop antennas, or saddle 
antennas. 

[0021] 

A schematic representation of some of the components of an NMR logging tool 30 
is illustrated in FIG. 2, which shows a first centralized magnet or magnet array 36 and 
an RF antenna J7, which may be a suitably oriented coil or coils. FIG. 2 also illustrates 
a general representation of closely-spaced cylindrical thin shells, 38-1 , 38-2 . . . 38- 
A/, that can be frequency selected in a multi-frequency logging operation. These thin 
shells are resonance regions where magnet 36 has a radial gradient in the field 
amplitude. One such device is disclosed in U.S. Patent No. 4,710,713 issued to 
Strikman. In FIG. 2, another magnet or magnet array 39 \s shown. Magnet array 39 
may be used to pre-polarize the earth formation ahead of the investigation region as 
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the logging device 30 is raised in the borehole in the direction of arrow Z. Examples of 





NMR logging tools having pre-polarization magnet arrays are disclosed in U.S. Patent 
Nos. 5,055,788 and 3,597,681 . 



pulses and for receiving and storing the spin echoes. One skilled in the art would 
appreciate that any other suitable circuitry could be used without departing from the 
scope of the invention. In FIG. 3, a downhole processor 210 has associated memory, 
timing, interfaces, and peripherals (not separately shown), as known in the art. The 
processor subsystem 210 is coupled with telemetry circuitry 205 , for communication 
with a processor (not shown) at the earth's surface. The pulse forming circuitry 
includes a variable frequency oscillator 220 which, under control of processor 210 , 
produces radio frequency (RF) signals at the desired frequencies. The output of 
oscillator 220 is coupled to a phase shifter 222 , which permits control of pulse 
phases, and then to a modulator 230 , both of which are under the control of 
processor subsystem 210 . The phase shifter Z?2and modulator 230 can be 
controlled, in a manner known in the art, to produce the desired pulse phases of RF 
field. The output of modulator 230 is coupled, via a power amplifier 235 , to the RF 
antenna 240 . A Q-switch 250 can be provided to damp the RF antenna system to 
reduce antenna ringing. The antenna 240 is also coupled with a receiver section via 
duplexer 265 , the output of which is coupled to receiver amplifier 270 . The duplexer 
265 protects the receiver amplifier 270 from the high power pulses which pass to the 
RF antenna 240 during the transmitting and damping modes. During the receiving 
mode, the duplexer 265 acts as a low impedance connection from antenna 240 to the 
receiver amplifier 270 . The output of receiver amplifier 270 is coupled to a dual 
phase-sensitive detector 275 , which also receives, as a reference, a signal derived 
from the oscillator signal. The detected output is coupled to analog-to-digital 
converter 280, the output of which is a digital version of the received nuclear 
magnetic resonance signal. Although the logging device or tool 30 \s shown as a 
single body in FIG. 1, it may alternatively comprise separate components, and the tool 
may be combinable with other logging tools. Also, while a wireline is illustrated, 
alternative forms of physical support and communicating link can be used, for 
example in a measurement while drilling system. 



[0022] 



FIG. 3 illustrates a schematic of a circuitry of an NMR tool for producing the RF 
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[0023] Several NMR parameters may be measured that can be used to derive formation 
properties. Most NMR logging operations measure the spin-lattice (longitudinal) 
relaxation times ( 7" ) and/or spin-spin (transverse) relaxation times ( T ^ ) of 
hydrogen nuclei. In addition, some NMR logging tools may provide a ratio of T I T 
directly, and other NMR tools may provide diffusion constants ( D). 

[0024] Various pulse sequences are known in the art for measuring the NMR relaxation 
times. For example, T (longitudinal) relaxation may be measured using an 
inversion-recovery (IR), saturation recovery (SR), or steady-state approach, or any of 
their derivatives. T ' relaxation is often measured from a train of spin-echoes that are 
generated with a series of pulses such as the Carr-Purcell-Meiboom-Cill (CPMG) pulse 
sequence or some variant of this. The CPMG pulse sequence is well known in the art. 
See Meiboom, S., Gill, D., 1958, "Modified Spin Echo Method for Measuring Nuclear 
Relaxation Times," Review of Scientific Instruments, 29, 688-91. CPMG pulse 
sequence generates a train of spin echoes, whose amplitudes exponentially decay as a 
function of time. The exponential decay life time is referred to as a transverse (spin- 
spin) relaxation time, T ' . Thus, T ^ measurements are accomplished by analyzing 
the amplitudes of spin echoes thus obtained. 

[0025] 

NMR measurements of diffusion constants are accomplished in the presence of 

magnetic field gradients. Magnetic field gradients produce different strengths of 

magnetic field at different locations. The different magnetic field strengths manifest 

themselves as different Larmor frequencies in the detected signals (because w = y B 

o 

). If a delay time is inserted between the pulse and the start of signal collection, 

o 

nuclear spins can diffuse from location A to location B, i.e., from one Larmor 
frequency to another, with a concomitant loss of phase coherence. This results in a 
decrease in the magnitudes of the detected signals. The diffusion induced decay of 
the NMR signal depends on the square of the echo spacing and on the magnetic field 
gradient, because the decay rate is governed by: 

where D is the diffusion constant, G is the gradient of the static magnetic field, y is 

the gyromagnetic ratio, and t is the echo spacing. See, e.g ., C.P. Slichter, Principles 

e 

rd 

of Magnetic Resonance, 3 Edition, Springer Verlag, Berlin, Heidelberg, New York, 
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1990, chapter 8 and appendix G. 

[0026] Once NMR data, which include information on T , relaxation times, T I 

ratio, or diffusion rates D , are collected, they are analyzed with an inversion method 
to derive a model of the earth formations. For the derived model, properties of the 
formation then can be calculated. Many methods are known in the art for inverting 
NMR data, including the Window Processing (WP) method disclosed in U.S. Patent No. 
5,291 ,1 37 issued to Freedman. However, when different fluids have similar NMR 
parameters, for example, same or overlapping T ^ distributions, most inversion 
methods cannot separately characterize these fluids. 

[0027] Freedman et al. disclosed a magnetic resonance fluid (MRF) characterization 

method that is capable of distinguishing different fluids even if they have overlapping 
NMR parameters (e.g., overlapping T ' distributions). See Freedman et al., " A New 
Method of Fluid Characterization in Reservoir Rocks: Experimental Confirmation and 
Simulation Results , "Paper SPE 6321 4 presented at the 2000 Society of Petroleum 
Engineers Annual Technical Conference and Exhibition, and Freedman et al., " Field 
Applications of a New Nuclear Magnetic Resonance Fluid Characterization Method " 
paper SPE 71713 presented at the 2001 Society of Petroleum Engineers Annual 
Technical Conference and Exhibition. In addition, U.S. Patent 6,229,308 Bl issued to 
Freedman discloses the application of this method in formation evaluation. 

[0028] The MRF method is based on two key concepts: (1 ) a new microscopic CVM 
(constituent viscosity model) that relates NMR relaxation times and molecular 
diffusion coefficients in crude oils, and (2) a new multi-fluid relaxation model. The 
MRF method provides a detailed formation evaluation of the near wellbore region 
investigated by modern NMR logging tools. The information provided by MRF includes 
flushed-zone fluid saturations and volumes, total and bound-fluid porosities, bulk 
volumes of hydrocarbon saturations, oil viscosities, and hydrocarbon-corrected 
permeabilities. 

[0029] CVM re | a tes individual constituent diffusion-free relaxation times and diffusion 
coefficients to a distribution of constituent viscosities. The constituent viscosities are 
molecular variables that are analogous to the "friction coefficients" used in Langevin 
equation models of Brownian motion in viscous media. Before looking at the 
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correlation between proton relaxation times and diffusion constants and the viscosity 
of a mixture, it is helpful to look at the simplest situation first, i.e., a pure liquid. 

[0030] In pure liquids, Bloembergen et al. in " Relaxation Effects in Nuclear Magnetic 

Resonance Absorption , "Physical Review (1948), 73, no. 7, pp. 679-712, showed that 
proton relaxation times (e.g., T ) depend on the temperature and the liquid viscosity 
according to: 

where a is a constant, T is the spin-spin (transverse) relaxation time, Tis the 

temperature of the liquid in degrees Kelvin, and n is the viscosity of the liquid. This 

o 

relationship is also valid for T , provided that the "fast motion" condition, w Q T < 
1 , is satisfied, where t is a rotational (i.e., molecular tumbling) correlation time and u) 
is the Larmor frequency. In pure liquids, the relationship between diffusion 

o 

constants and viscosities is described by the Einstein-Stokes model, which predicts 
that diffusivity is related to temperature and viscosity according to the equation: D = 
kT/6 n rj R , where k is the Boltzmann constant, R is the radius of the spherical 
particle, n is the viscosity of the fluid, and 7"is the temperature in degrees Kelvin. This 
relationship may be simplified as D = bT/ n . where b '= k/6 tt R. 



[0031] 



The dynamic properties of each individual constituent in a mixture are similar to 

those of a pure liquid. Thus, a similar relationship exists between the relaxation time 

( T , ) and the diffusion constant ( D , ) of the k-th constituent in a mixture of 
2o,k o,k 

oils according to CVM. The subscripts " o,k " denotes the parameter is for the k-th 
constituent in the oils. The key CVM equations for a hydrocarbon mixture with n 
components are: 

^^rf 1 -fo)L (3) 

where T , and D , are the relaxation time and the diffusion coefficient, 
2o,k o,k 

respectively, for the k-th molecular constituent in a hydrocarbon mixture (e.g., a 
crude oil), T\s the sample temperature in degrees Kelvin, f(COR) , as will be explained 
below, is an empirically determined "universal function" of gas/oil ratio (GOR) and is > 
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[0033] 




1 , a and b are constants, rj is the constituent viscosity for the k-th molecular 

K 

constituent, f, is the proton fraction of each of the mixture constituents, and rj is 
k o 

the macroscopic viscosity ( rj ) of the mixture. Equation (3) indicates that the 

o 

macroscopic viscosity of the mixture equals the logarithmic mean (or geometric mean) 

of the constituent viscosity distribution (( rj ).._). 

k LM 

[0032] The factor f(GOR) is included in equation (1) because it has been shown that GOR 
is an important parameter in determining the relaxation time dependence on viscosity 
and temperature. See Lo et al., Relaxation Time And Diffusion Measurements of 
Methane And N-Decane Mixtures , The Log Analyst, pp. 43-46, November-December, 
1 998; see also U.S. Patent No. 6,229,308 B1 . Note that the f(GOR) function is always 
equal or greater than 1 . For dead oils (which are oils having no solution gas) and 
mixtures that do not contain any solution gas, f(GOR) = f(0) = 1 . (Freedman, R M et al., 
" A New Method of Fluid Characterization in Reservoir Rocks: Experimental 
Confirmation and Simulation Results /'paper SPE 6321 4 presented at the 2000 SPE 
Annual Technical Conference and Exhibition, Dallas, 1-4 October; hereafter referred 
to as "Freedman et al."). The constant a in equation (1) has been empirically 

determined to be about 0.004 s • cp • K 1 and this value has been shown to be 
applicable to substantially every constituent in crude oils. See Freedman et al. The 

-8 

constant b \r\ equation (2) has been empirically determined to be about 5.05 x 10 

cm 2 • s 1 • cp • K V This value for constant b is also universally applicable to 
constituents in both dead and live hydrocarbon mixtures and crude oil. See Freedman 



et al. 



Combining equations (1) and (2) provides the expression: 

Equation (4) indicates that there is a linear relationship between the diffusion constant 

( D , ) and relaxation time ( T , ) for each of the constituents in the mixture. This 
o,k 2o,k 

linear relationship is similar to the relationship described above for the pure liquids. 

Furthermore, CVM should converge with a pure liquid model when there is only one 

constituent in the mixture (i.e., pure liquid, n = k =1). Indeed, when n = 1 , 

o o 

equations (1) and (2) reduce to a form identical to that for pure liquids as described by 
Bloembergen et al. 
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7-/ (MB) J? 0 



and Stokes Einstein model 



[0034] The validity of equations (1) - (4) has been shown (as disclosed in Freedman et al.) 
using experimental measurements of T and £> distributions on live and dead 
hydrocarbon mixtures and crude oils. Equations (1) and (2) provide some interesting 
insights into the nature of the molecular interactions in hydrocarbon mixtures such as 
crude oils. First, each molecular constituent retains its separate properties in the 
mixture. Second, the relationships between NMR relaxation times, molecular diffusion 
coefficients, and constituent viscosities for each constituent in a mixture are the same 
as that for a pure liquid, except that the constituent viscosity rather than the pure 
fluid viscosity (macro viscosity) determines the relaxation time and the diffusion 
coefficient of each constituent. Third, all unknown factors that might affect molecular 
interactions, e.g., molecular composition and molecular sizes, are contained in the 
constituent viscosity terms in this model, and there is no need to determine these 
factors. In addition, equation (4) predicts that the T ^ and D distributions in 
hydrocarbon mixtures and crude oils are not independent. 

[0035] The importance of the CVM for fluid characterization can be appreciated from 

equation (4), which correlates diffusion-free relaxation times with molecular diffusion 
in crude oils. This link reduces the number of unknown parameters in the NMR 
multifluid relaxation model and results in robust and accurate recovery of oil T ' and 
D distributions using the MRF method, even when the oil and connate water (or brine) 
T ^ distributions are broad and largely overlap one another. 

[0036] The MRF method uses a general relaxation model for a formation containing brine, 
hydrocarbons, and oil-based mud filtrate (OBMF). Mud having hydrocarbons as the 
continuous phase is commonly used in drilling the borehole and it invades porous 
formations as OBMF. In the absence of magnetic field gradient or diffusion, the 
amplitudes of the spin echoes in NMR measurements decay exponentially as a 
function of time ( 7" or T ' ). In the presence of a magnetic field gradient, molecular 
diffusion also contributes to the reduction of the measured NMR signal amplitudes. 
Let 
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represent the amplitude of the j-th echo acquired during measurement p. Consider 
the following general relaxation model, 



•(-£)) 



1-exp 



(5) 

where the first, second, and third terms are connate water, native hydrocarbon, and 
OBMF signals, respectively. This three-phase model explicitly accounts for 
contributions from individual constituents in the connate water and hydrocarbon 
phases, but only assumes an average relaxation time distribution in the OBMF phase. 
The simplified term for OBMF is justified because experimental measurements in 
OBMF have shown that NMR relaxation time distributions for individual constituents in 
typical OBMF are very narrow and can be described by a single decay exponent. If such 
simplification is not desirable, one may include a summation function, like that for 
connate water or native hydrocarbon, for the OBMF. In addition, one skilled in the art 
would appreciate that this general relaxation model may be modified, for example, to 
remove the OBMF term, if justified (e.g., investigation of core samples removed from 
formations, or LWD-NMR logs in which mud infiltration is negligible), or to include gas 
in the fluids as disclosed in Appendix A of U.S. Patent No. 6,229,308 Bl. 

[0037] The apparent transverse (spin-spin) relaxations ( T ^ ) in any of the three phases 

modeled in equation (5) include actual transverse relaxation and signal decay due to 

self diffusion of molecules in the static magnetic field gradient G . For unrestricted 

P 

diffusion in a uniform magnetic field gradient in the connate water, the apparent 
transverse relaxation rates can be written as, 

Here, T in the first term on the right-hand side are a set of relaxation times that 
2 1 1 

represent the sum of surface and bulk relaxation of the connate water phase. The 
second term is the contribution to the relaxation rate from diffusion, where 

is the proton gyromagnetic ratio and D ^ (T) is the temperature dependent self 
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diffusion coefficient of water in units of cm ^ /s. Note that equation (6) assumes an 

unrestricted diffusion and a uniform magnetic field gradient C . One skilled in the 

p 

art will appreciate that corrections to D CO for the effects of restricted diffusion and 

w 

to G for the effects of internal rock gradients can be applied if appropriate, such as 
P 

when the formations have minerals having high magnetic susceptibility. 
[0038] Similarly, the apparent transverse relaxation rates in the native oil 

can be written in the form, 

5 ° M « 

where T ( n ) is the bulk relaxation time associated with amplitude b in the 
2,o k k 

hydrocarbon relaxation time distribution, and D (rj ) \s a viscosity dependent 

o k 

diffusion constant. 

The measured viscosity ( q ) reflects a macroscopic transport property of the 

o 

crude oil that determines its flow properties and is the quantity that is used in 

hydrodynamic transport equations like the Navier-Stokes equation. Morriss etal. 

showed that, for a suite of dead (i.e., not containing dissolved solution gas) crude oils, 

there exists a strong correlation between the logarithmic mean relaxation times of 

their constituents and the measured viscosities. See Morriss et al., " Hydrocarbon 

Saturation And Viscosity Estimation From NMR Logging in The Be/ridge Diatomite , 

"Paper C presented at the 35th Annual Meeting Of The Society Of Professional Well 

Logging Analysis, 1 994. The macroscopic viscosity ( n ) of live crude oils is 

o 

empirically related to the logarithmic mean 

of the transverse relaxation time distributions by an equation of the form, 



[0039] 



1 _ anJiOOK) 



a cv» (8) 



where a is an empirically determined constituent constant that has been determined 
by Looyestijn (" Determination of Oil Saturation From Diffusion NMR Logs ," Paper SS 
th 

presented at the 37 Annual meeting of the Society of Professional Well Log 
Analysts, 1 996) to be around 250 (i.e., a s 250 Ks cp ), for 

in seconds and rj in centipoise and Tthe temperature in degrees Kelvin. Thus, c = 
o 
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a I 7"in dead oil (no dissolved gas) or 

C " T 

for live oil. The empirically derived function f(COR) accounts for live oils (those 
containing dissolved solution gas) and has been discussed by Freedman et al . in 
paper number 63214 entitled "A New NMR Method of Fluid Characterization in 
Reservoir Rocks: Experimental Confirmation and Simulation Results" presented at the 
2000 Society of Petroleum Engineers Annual Technical Conference and Exhibition 
meeting. 

[0040] The q . terms in equation (7) are microscopic viscosities that reflect the complex 
composition of crude oils. Analogously with equation (8), the constituent viscosities 
are assumed to be related to the components in the relaxation time distribution via 
the same equation, 

1 *TlJ<GOR) 

(9) 

[0041] The dependence of the relaxation times on viscosity and temperature in equations 
(8) and (9) is consistent with the experimental observations and theoretical 
predictions of Bloembergen, Purcell, and Pound, Relaxation Effects in Nuclear 
Magnetic Resonance Absorption , Physical Review, vol. 73, no. 7, pp. 679-712, 1948. 

[0042] Similar to the Stokes-Einstein equation, the self-diffusion constants for the crude 

oils, D , are assumed to have the same dependence on Tf r\ , . Therefore, for crude 
o k 



oils, 



7* 



(10) 

where b is a constant, D is the measured crude oil diffusion constant in cm 2 /s, and 
7"is the temperature in degrees Kelvin. The empirical constant for crude oils, 

5.05* lO" 3 cot VW' 1 

, is given by Freedman etaL in paper number 63214 entitled "A New NMR Method of 
Fluid Characterization in Reservoir Rocks: Experimental Confirmation and Simulation 
Results, "presented at the 2000 Society of Petroleum Engineers Annual Technical 
Conference and Exhibition meeting. 

[0043] Analogously to the relationship between the macroscopic diffusion constant ( D ) 

o 
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and the macroscopic viscosity ( rj ), the microscopic constituent diffusion constants 

o 

( D ( n )) of crude oil mixtures are related to the microscopic constituent viscosities 
o k 

( rj ) according to the following equation: 

IN 

DM -- *H*.0-> . 

GO 

Equation (1 1 ) implies that there exists at the molecular level a distribution of diffusion 
constants in the crude oil mixture analogous to the distribution of relaxation times. 
This is the key concept of the CVM (constituent viscosity model). These distributions 
of diffusion constants and relaxation times can be derived from the measured NMR 
data by iteratively fitting a multifluid model as shown in equation (5) to these data 
using the method disclosed in U.S. Patent No. 6,229,308 Bl issued to Freedman. 

[0044] The MRF method inverts suites of NMR data using a multifluid relaxation model as 
shown in equation (5). In this model, different fluids that have overlapping T 
distributions are distinguished based on differing diffusion constants. This is feasible 
because water molecules typically diffuse much faster than oil molecules, and much 
slower than gas molecules. Consequently, NMR diffusion measurements can provide a 
means for detection and differentiation of reservoir fluids, even if other NMR 
parameters cannot. 

[0045] jo use t h e contributions of differing molecular diffusion rates to resolve the fluid 
composition where there are overlapping T ^ distributions, the MRF method uses a 
suite of spin-echo measurements which have been acquired with different echo 
spacings, polarization times, applied magnetic field gradients, and/or numbers of 
echoes. The suite of NMR data may be acquired with a conventional CPMG-pulse 
sequence or a new type of "diffusion-editingTDE") spin-echo sequence tailored for 
fluid typing. See M.D. H u rlimann et al., " Diffusion-Editing: New NMR Measurement 
of Saturation and Pore Geometry " paper presented at the 2002 Annual Meeting of the 
Society of Professional Well Log Analysts, Osio, Japan, June 2-5; see also , U.S. 
Application Serial No. 09/723,803 filed on November 28, 2000 by H u rlimann, 
entitled " Simpler and More Robust Hydrocarbon Typing with Pulsed NMR This 
application is assigned to the assignee of the present invention and is hereby 
incorporated by reference. DE sequences are similar to the CPMG sequences except 
that the initial two echoes are acquired with longer echo spacings and the third and 
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[0048] 




subsequent echoes are acquired with shorter echo spacings. Diffusion information is 
encoded during the acquisition of the first two echoes, whereas the third and 
subsequent echoes can provide bulk and surface relaxation time information with 
little, if any, attenuation of the signal by diffusion. 

[0046] The suite of spin-echo measurements are then used in inversion analysis to fit the 
general multifluid relaxation model as shown in equation (5). Any inversion method 
known in the art may be used, including the Window Processing (WP) method 
disclosed in U.S. Patent No. 5,291 ,137 issued to Freedman. This patent is assigned to 
the assignee of the present invention and is hereby incorporated by reference. For 
more detailed discussion of the MRF methods see Freedman et al., " A New Method of 
fluid Characterization in Reservoir Rocks: Experimental Confirmation and Simulation 
Results ," paper SPE 632 1 4 presented at the 2000 SPE Annual Technical Conference 
and Exhibition, and Freedman et al., " Field Applications of a New Nuclear Magnetic 
Resonance Fluid Characterization Method " paper SPE 71713 presented at the 2001 
SPE Annual Technical Conference and Exhibition. 

[0047] As stated above, for the diffusion constants to be measurable, the nuclear spin 
relaxing processes (longitudinal or transverse relaxation) should be slower than the 
diffusion process. Otherwise, the NMR signals would have diminished too much before 
any appreciable diffusion has occurred. In other words, if the NMR signal decays too 
quickly such that the decay becomes comparable with or faster than the signal decay 
due to the diffusion process, the diffusion constants can no longer be determined. 
Without the diffusion constants, the fluids with overlapping T ' distributions cannot 



be separated by the conventional MRF method. 



The lower limits of T or , above which the MRF method is useful, are even 

TM 

higher for logging-while-drilling (LWD-NMR ) tools because these tools typically 
have low magnetic field gradients. As is evident from the above discussion, the 

magnitudes of the diffusion-induced relaxation depend on the strengths of the 

TM 

magnetic field gradients. With low-gradient NMR tools (e.g., LWD-NMR or MRWD 
TM 

), it requires longer delay times (e.g., the inter-echo delays) for the diffusion 
information to become detectable in the measurements. The required longer delay 

time makes the diffusion measurement only possible for fluids with relatively long 7"^ 
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or T times. For example, a typical LWD-NMR tool operating at a gradient of 4 G/cm 
(Gauss per centimeter) would need about 6 times longer delays, as compared with a 
TM 

typical MRX tool operating at a gradient of 25 G/cm, for the molecular diffusion to 
be measurable. Thus, using a low-gradient NMR tool, the diffusion-based MRF 

method is useful only for oils with long 7" or T ^ values, i.e., free fluids, but not for 

heavy oils. 

[0049] Embodiments of the present invention are capable of extending the applicability 
of the MRF method to situations where short NMR relaxation times do not permit 
accurate determination of the diffusion constants. Embodiments of the invention use 
other characteristics of the fluids, for example, different T ^ I T ^ ratios ( \ in 
equation (5)), to differentiate fluids having similar or overlapping T ^ distributions. T 
^ I T ' ratios are readily measurable using typical NMR logging tools. T ^ / T ^ 
contrasts (different 7" / T ' ratios between different fluids) have been used in the art 
for other purposes, such as "polarization correction," which is important for 
measurements acquired using pulse sequences that have short wait times (WT). See 
U.S. Patent No. 6,229,308 Bl , issued to Freedman. 

[0050] In free fluids, T ' equals T ^ except for very viscous fluids (e.g., tar) and/or for 

very high Larmor frequencies (usually not encountered in logging tools). Thus, in free 
fluids, the 7" / T ratio equals 1 . The porous spaces in the formations often are 
water wet. As a result, oils (hydrocarbons) are repelled from the surfaces (walls) of the 
pores, and they behave like free fluids without boundary effects. In contrast, water 
interacts with the surfaces of the rock grains. Consequently, each water molecule 
would spend some time as a surface-bound water molecule and other time as a free- 
water molecule. Surface-bound water molecules tend to have faster transverse (spin- 
spin) relaxation, i.e., short T ' , while free-water molecules would have longer T ^ 
values that approach the 7" values. Therefore, on average, water molecules in the 
formations have shorter (average) T ^ values than those of free fluids, e.g., oils 
(hydrocarbons). Consequently, for fluids in the formations, there is a T I T ' 
contrast between oil and water: T I T (oil) 1, and T ^ / T ^ (water) > 1, e.g., 1.5 



[0051] 



3.0 in sandstones. 



A T ^ I 7" 2 contrast between oil and water of 1 .5 does not seem large. However, 
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this ratio appears as an exponent in an exponential function that describes the echo 
train amplitudes ( see e.g., equation (5)). Consequently, 7" / T ^ contrast of 1 .5 
provides an easily detectable amplitude difference between two fluids with identical T 
^ times. 

[0052] The T I T contrast can be exploited using several NMR methods, for example, 
saturation recovery (SR), inversion recovery (IR), or steady state method in 
combination with conventional CPMG pulse sequences. Among various methods for T 
i evaluation, IR (which monitors the growth of the signal magnitude from 1 to +1) has 
more-dynamic range for signal detection than SR does (which monitors the growth of 
the signal magnitude from 0 to +1) (see FIG. 4A). In FIG. 4A, curve (a) depicts the 
exponential growth of an NMR signal magnitude after the saturation power has been 
turned off, and curve (b) depicts the exponential growth of an NMR signal magnitude 
after a 1 80-degree inversion pulse. It is apparent from FIG. 4A that curve (b) would 
provide at least two fold better sensitivity in monitoring the exponential growth than 
would curve (a). As a practical matter, the measurements cannot start immediately 
after a saturation pulse is turned off. A delay time is required for instrument circuitry 
to equilibrate before the detection can begin. This delay further degrades the 
sensitivity of the SR approach because most dramatic magnitude changes are lost 
during the delay time. For this reason, IR typically produces data that are about 3 fold 
better quality than those from SR. Therefore, the following discussion will use IR as an 
example. However, one skilled in the art would appreciate that the following 
discussion applies equally well to all other methods that can produce 7" / T ^ ratios, 
including SR. 

[0053] Using inversion recovery (IR) as an example, a pulse sequence for measuring the T 

I T ratio would be: (WT 1 80 RT CPMG) , , as illustrated in FIG. 5A, where WT is 
12 n 

the wait time (also called "polarization time") for spins to relax back to the steady 
state, RT is the recovery delay time during which the signals exponentially grow back 
to the steady state magnitudes, and CPMG is the typical CPMG pulse train (i.e., a series 
of 1 80-degree pulses with a fixed delay time, TE, between them). This pulse sequence 
will be referred to as an "IR-CPMG" sequence. 

[0054] 

In FIG. 5A, the first 1 80-degree pulse inverts the magnetization (i.e., onto the Z 
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axis). After /?rdelay the signal may be detected with a 90-degree pulse. As shown in 
FIG. 5A, the signals are actually detected with a CPMG sequence as a train of spin 
echoes. The detected spin echo magnitudes will be influenced not only by the T ' 
relaxation process, but also the T relaxation process, which is encoded by the 
various /?7~delay times. The IR-CPMG sequence as shown in FIG. 5A would detect 
signal components that have relaxed back to. the Z axis after /?rdelay. Herein, the "IR- 
CPMG" would be used to refer to the sequence shown in FIG. 5A and any variant 
sequence which substantially include this sequence and achieve the same purpose. 
FIG. 5B shows a CPMG sequence coupled with a saturation recovery (SR) sequence, 
instead of an inversion recovery (IR) sequence. This sequence would be referred to as 
an "SR-CPMG" sequence. The SR-CPMG sequence as shown in FIG. 5 B would detect 
signal components that have relaxed back to the Z axis after tfrdelay. As stated 
above, the SR-CPMG may be used instead of the IR-CPMG in the embodiments of the 
invention. 

[0055] Ideally, WT in any pulse sequence should be long enough for all spins to relax 
back to the steady state (i.e., fully polarized by the static magnetic field) before the 
next pulse cycle. However, time economy often dictates that a shorter WT is used, 
which may result in incomplete magnetization being detected. If insufficient WT is 
used, polarization of the spins by the static magnetic field may not be complete. The 
detectable signal intensity under this condition is given by: 

where / denotes the detectable signal intensity when the spin is fully polarized by 
o 

the static magnetic field (or fully relaxed back to the steady state). 



[0056] 



The measured signal amplitude ( A . for the i-th echo from the CPMG train) using 
an IR-CPMG sequence would be: 



4 ^l-S* J«~"% (1 1) 



where A is the full signal amplitude after full repolarization along the static magnetic 
field, RTis the inversion recovery time in the IR-CPMG sequence (see FIG. 5A), 7"£is 
the delay time between the 1 80-degree pulses in the CPMG sequence, T and T ' are 
the longitudinal and transverse relaxation times, respectively, and 5 is defined as: 
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where IE\s the inversion efficiency. The inversion efficiency, /£, is 1 .0 under ideal 
condition. However, /£is around 0.8 for a typical well logging tool. 

[0057] The above formulation assumes that signal decay due to diffusion is negligible, 
and the measured spin echoes decay with true T . In addition to the T ^ decay 
during the CPMC pulse train, the measured spin echo magnitudes are also influenced 
by the /?7~delay times in the IR-CPMC sequence ( 7" relaxation). Alternatively, if RT is 
kept constant, the signal magnitudes would be influenced by T relaxation. 

[0058] Referring to curve (b) in FIG. 4A, it is apparent that there exists a critical RT, for 

which the recovering signal would vanish and become undetectable (i.e., A . = 0). This 

/?roccurs at where curve (b) intercepts the time axis in FIG. 4A. Because different 

fluids have different T values, each fluid may have a different flrtime when the 

signal would become undetectable. Referring now to FIG. 4B, if a fixed RT delay time, 

e.g., T , is used in the IR-CPMG sequence, some signals would have recovered 
1 c 

sufficiently to produce a positive signal magnitude (curve (c)), some might become 
undetectable (curve (d)), and some may produce a negative magnitude (curve (e)). 

[0059] According to equation (11), A . = 0 when T = RT I ln( 5). Thus, T^ = RT I 

0.59 for 5= 0.8, or = RT I 0.69 for 5= 1.0. Signal components with T ^ < T^ 

(curve (c) in FIG. 4B) will appear with A > 0, and signal components with T > T 

(curve (e) in FIG. 4B) will appear with A < 0, while signal components with T = T 

i 1 lc 

(curve (d) in FIG. 4B) will be undetectable. 

[0060] FIG. 4B illustrates the node where A . = 0 in T measurement. This node can also 
be located using the first echoes in the CPMG measurements using the IR-CPMG 
sequence shown in FIG. 5A having different /?rdelay times. In well logging, a train of 
CPMG echoes are typically recorded for each flrdelay. The node in the measured data 
where the intensity is zero will depend not only on T but also on T . Therefore, the 
position of this node (where A . = 0) in the T ^ spectrum is determined by the T ^ / T 
^ ratio. By recording echo trains for several RT, it is possible to invert the data for two 
components with two T ^ / T ^ ratios in the neighborhood of this node. The simplest 
measurement suite would consist of two measurements: one with a standard burst 
CPMG without the inversion recovery sequence and the other with an IR-CPMG with RT 
< WT. 
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[0061] For data inversion, if the T / T ratio for water is known, inversion routines, 
such as the Window Processing method as disclosed in U.S. Patent No. 5,291 ,1 37 
issued to Freedman, for normal CPMC data suites may be used. Only the data model 
for the experimental data has been changed to that described above. Data inversion 
would yield T ^ distributions for oil and water directly. From these T ^ distributions, 
porosity fractions (water-filled porosity, oil-filled porosity, and total NMR porosity), 
volume fractions (bound fluid volume and free fluid volume), constituent oil 
viscosities, viscosity or diffusion constants of the oil mixture or brine, and 
permeabilities can then be calculated. Methods for deriving formation properties from 
the T ' distribution are known in the art. For example, U.S. Patent 6,229,308 B1 
provides details on how this might be achieved. 

[0062] FIG. 6 illustrates key elements in the embodiments of the present invention. First, 
a suite of (i.e., a plurality of sets of) NMR measurements are acquired (shown at 81 ) 
with at least one parameter selected from wait time or polarization time (WT), recovery 
time (RT), inter-echo delay time (TE), and number of echoes (NE) varied between 
different sets of measurements, wherein each set of measurements includes 
information on the longitudinal ( T ) and transverse ( T ^ ) relaxation process. For 
example, the measurements may be acquired with the IR-CPMG or SR-CPMG sequence 
(as shown in FIGs. 5A and 5B), or any variation of these sequences so long as the 
collected data include information on both T and T . The parameter is preferably 
so varied that the T and/or T ' process affects each NMR data set to a different 
extent. The sets of NMR data may be acquired with any NMR logging tool, wireline or 
logging-while-drilling. In addition, these measurements may be made with core 
samples removed from the formations in the laboratory. 

[0063] 0nce NMR data are ^11^^ a formation model is generated (shown at 82) for 
inversion analysis. This model may include oils, brine, OBMF, and/or gas. In addition, 
within each fluid type (e.g., oils), the model may include multiple components to 
account for the constituents in the mixture. The model may include estimated signal 
magnitude, longitudinal relaxation time, transverse relaxation time, ratio of 
longitudinal to transverse relaxation time, diffusion constant, and viscosity for each 
constituent. This model is then optimized (shown at 83) to fit the collected data. The 
optimization process typically includes calculating a synthesized NMR data set based 
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on the model; comparing the synthesized NMR data set with the NMR measurements; 
and adjusting the model and repeating the calculating and comparing processes, if 
the difference between the synthesized NMR data set and the NMR measurements is 
greater than "a minimum." One skilled in the art would appreciate that "a minimum" 
for optimizing a formation model as used herein may be a true minimum of the 
difference between the calculated values and the experimental values. Alternatively, "a 
minimum" for the purpose of formation optimization may be a preset value, below 
which the difference between the model and the formations is considered acceptable. 
The optimization process may use any algorithm known in the art, for example 
PowelTs nonlinear Constrained Minimization Algorithm cited in U.S. Patent No. 
6,229,308 Bl , or Monte Carlo simulation. 

[0064] The optimized formation model is then output (shown at 84) with the computed 
parameters for each constituent fluid. These parameters may include 7" or T 
relaxation time, diffusion constants, and constituent viscosity. From these constituent 
parameters, the desired formation properties can then be calculated (process 85). As 
stated above, the formation properties may include: porosity fractions (water-filled 
porosity, oil-dilled porosity, and total NMR porosity), volume fractions (bound fluid 
volume and free fluid volume), constituent oil viscosities, viscosity or diffusion 
constants of the oil mixture or brine, and permeabilities. 

[0065] The utility of embodiments of the present invention will now be illustrated with 

the following examples, in which computations are performed with a software 

TM 

package such as that sold under the trade name of MatLab by The MathWorks, 



First, two identical T distributions for oil and water (representing the most 

difficult case to resolve) were created with T values in the range from 2.5 to 25 ms. 

The oil filled porosity is 1 0 pu, the water-filled porosity is 20 pu. The T I T ^ ratio 

was set to 1 .0 for the oil component and 1 .5 for the water component. These 

distributions of T were observed with 3 suites of 5 IR-CPMG or SR-CPMG and 1 
2 

CPMG sequences, having the following parameters: 



Natick, MA. 



[0066] 
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# * 







Suite 1 : Long WT IR Suite 






* IE 


0.8 


0.8 


.0.8 


0.8 


0.8 


NA* 


WT 


10,000 


10,000 


10,000 


10,000 


10,000 


10,000 


RT 


1 


Jo 


10 




100 


NA* 


TE 


0.5 


0.5 


0.5 


0.5 


0.5 


0.5 


NE 


500 


500 


500 


500 


500 


500 



Parameter is not applicable, since this represents the normal CPMG sequence. 







Suite 2: Long WT SR Suite 






IE 


0.0 


0.0 


0.0 


0.0 


0.0 


NA 


WT 


10,000 


10,000 


10,000 


10,000 


10,000 


10,000 


RT 


1 




10 




100 


NA 


TE 


0.5 


0.5 


0.5 


0.5 


0.5 


0.5 


NE 


500 


500 


500 


500 


500 


500 








Suite 3: Optimal WT IR Suite 






IE 


0.8 


0.8 


0.8 


0.8 


0.8 


NA 


WT 


10 




100 . 


jTooa -vjio.ooo 


10,000 


RT 


1 




10 




100 


NA 


TE 


0.5 


0.5 


0.5 


0.5 


0.5 


0.5 


NE 


500 


500 


500 


500 


500 


500 



Note that RT = 0 is used as a flag to describe a standard CPMG preceded by WT. All 
times are in milliseconds (ms). Each sequence was repeated only once with 500 
echoes collected (NE = 500). 

[0067] 

One pu of Gaussian noise was added to the experimental data. Then, 25 noise 

realizations were run for each suite. The data were then inverted using the nnlsQ 
TM 

function in MatLab . This function provides an inversion with positivity constraint. 
The data were either inverted directly or with window sums. The window sums 

contained the following echoes: 1 , 2-3, 4-6, 7-1 0, 11 -1 6, 1 7-25, 26-39, 40-60, 61 - 

91, 92-138, 139-208, 209-314, 31 5-500. The following tables show the results from 

the inversions (using "director "window sums" (WS) approach): 
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T 2 (ms) 





Oil 


Water 




1^=6.89 




Suite 


T 2]a (direct) T^WS) 


(direct T^(WS) 


1 


7 .3 ±2 .2 


7 .0 ±13 


2 


73 ±2.4 


6.9 + 1.2 


3 


6 .5 ±1.6 


12 ±L0 



Porosity 






Oil 


Water 




<Ho.io 




0.20 


Suite 


4> (direct) 




4> (direct) 


4>(WS) 


1 


0.11 ±0.02 


0.11 ±0.01 


0.20 ±0.02 


0.20 ±0.02 


2 


0.11 ±0.03 


0.11 ±0.03 


0.20±a03 


0.19±0.03 


3 


0.10 ±0.02 


0.10 ±0.01 


0.20 ±0.02 


0.20±0.02 



"WS"indicates inversions performed with window sums. 

[0068] These results show that the IR-CPMG suites (suites 1 and 3) provide better quality 
data than does the SR-CPMG suite (suite 2), as evidenced by the smaller standard 
deviations. As stated earlier, inversion recovery (IR) experiments monitor the growth 
of signals from 1 to +1 , while the saturation recovery (SR) experiments monitor the 
growth of signals from 0 to + 1 (see FIG. 6A). Thus, it is expected that the IR-CPMG 
suites would provide better quality data than would the SR-CPMG suite. 

[0069] A comparison between results from suites 1 and 3 reveals that suite 3, which has 
optimized WT values, produces better results. This is probably due to the spread of 
WT values that provide more accurate fitting of data to the model. This is so because 
the amplitude of each echo includes the term 

S-=l+jffi*|l-e j, 

which can be better fitted with WT values spanning the T value (i.e., from WT < T ^ 
to WT > T ), as compared with a single WT value in suite 1 . In addition, the 
parameters in suite 3 permit faster signal acquisition (than suites 1 and 2). The faster 
acquisition will allow more signal averaging/stacking within the same period of time, 
hence better signal-to-noise ratio (SNR). 
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[0070] The above simulation has assumed a T ^ / T ratio of 1 .0 for oils and a 7"^ / T ^ 

ratio of 1 .5 for water. To investigate the effects of 7" I T ' ratio mismatch on the 

results, simulations were run with various T ^ / T ' (water) ratios, ranging from 1.5 to 

2.4 in steps of 0.1. The same T ^ distributions for oil and water were used. For each 

ratio, 100 noise realizations were averaged. Results for porosities and T are 

21m 

shown in FICs. 7 A and 7B, respectively. 

[0071] As shown in FIG. 7A, the porosity results are relatively sensitive to the T I T ^ 

ratio. A change of 7" / T ' from 1 .5 to 2.4 resulted in an about 2-fold change in the 

porosities. In contrast, FIG. 7B shows that ^ 2 | m is ' ess sensitive to 7^ / T ^ changes. 

A change of T IT from 1 .5 to 2.4 resulted in less than 40% changes in the T , 
1 2 2lm 

values. 

[0072] Because the results depend very sensitively on the T I T ratio, this quantity 
needs either to be known a priori or it needs to be determined out of the data set as 
well. In principle, the data fit should be best for the correct 7" / T ratio, i.e., the 
error function (square sum of misfit between measurements and reconstruction) 
should be at a minimum for this value. In case of a linear model, this function can be 

written as: f= I ( mdata{ Mx spec)) ^ , where mdata are the measurements, M 
describes the linear model, spec is the vector of model parameters, and Mx spec 

describes the reconstructed NMR data based on the model. The summation goes over 

all measured data. 

[0073] FIG. 8 shows that there is, indeed, such a minimum that can be used to fix the T 
^ I 7" 2 ratio of water. The results in FIG. 8 are obtained with inversion of the window 
sums as described above. As can be seen from this inversion (which is averaged over 
1 00 noise realizations), the minimum is very flat and thus not well defined. However, 
such minimum may be better defined by averaging over large intervals with similar 
lithology (i.e., similar surface relaxivity). 

[0074] instead of using a common T ^ I T ^ ratio for each fluid averaged over a large 

depth interval as illustrated above, it is also possible to determine the 7" / T ratio 
for each depth level in the inversion using a non-linear fit routine, e.g., the inversion 
method of Freedman (G.B. Patent No. 2,338,068 or U.S. Patent No. 6,229, 308B1) or 
the above cited non-linear optimization routine by Powell. However the weak 



App_JD : 



10064914 



Page 25 of 41 




dependence of the IR-CPMC data suite on the T ^ I T ratio makes the accuracy of 
the determined value questionable. Thus, the inversion using the T ^ I T ratio for 
each depth level may not be more accurate than that using a common T / T ^ ratio. 

[0075] The above discussion presumes a homogeneous field. In an inhomogeneous field, 
the echo train decay constant may have contributions from both T ^ and T , and 
these contributions may change with the echo number. Thus, the effective T I T 
ratio for water might be echo number dependent for logging tools. In addition, it is 
assumed that the T / T ' ratio is constant for each fluid. For t values in the BFV 
(bound fluid volume) region, the 7" / T ratio is probably constant. However, this 
may not be true for all fluids. For example, for water in large vugs, surface relaxation 
will be unimportant and thus the ratio will approach 1.0. If this situation is suspected, 
it is advisable to use other types of measurements to confirm the results. In addition, 
if the T ' I T ratio changes with T ^ , the function should be known or be 
determined. This function could be determined from cores or from analysis of a 
similar non-reservoir interval (without oil). 

[0076] Advantages of the present invention include its applicability to spins with short T 

or T , where diffusion-based approach fails. In addition, with short T and T ^ , 
12 12 

the IR-CPMC sequence may be repeated quickly to gather good SNR data in a 
relatively short period of time. Furthermore, the invention is not limited by the types 
of NMR tools used for acquisition; it is applicable to a wireline NMR tool, a formation 
tester having an NMR module, or an LWD NMR tool. In addition, methods of the 
invention may be applied to NMR measurements obtained in a laboratory using a fluid 
sample, which may be any mixture of fluids whether formation fluids or not. 

[0077] The above discussion uses the T I T ' contrast as a parameter in the MRF 
method to separate fluids. Other parameters may be used in combination. -For 
example, diffusion contrast may be used together with the 7"^ / T ' contrast. These 
information may be included in one data set and simultaneously inverted for optimal 

SNRf<Jt 

for the answer products. If the diffusion is not taken into account, the echo spacing 
preferably should be short so that diffusion is negligible and the decay constant is the 
diffusion-free T value. 



AppJD=10064914 



Page 26 of 41 



[0078] While the invention has been described with respect to a limited number of 

embodiments, those skilled in the art, having benefit of this disclosure, will appreciate 
that other embodiments can be devised without departing from the scope of the 
invention as disclosed herein. Accordingly, the scope of the invention should be 
limited only by the attached claims. 



App_ID= 100649 14 



Page 27 of 41 



